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ABSTRACT 

A dynamic inventory control model, formulated in terms of 
discrete variable servomechanism theory, is examined. The 
behavior of the system is first analyzed by digital computer 
simulation. Both random and deterministic inputs are considered. 
Two double echelon models are then constructed and examined 
analytically utilizing transform theory. Finally, the proper- 
ties of forecasters of demand are analyzed for the cases where 
mean demand is constant and linear. Exponential smoothing 
techniques are compared with standard estimation procedures. 
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1. INTRODUCTION 



This paper will take for its starting point a dynamic 
inventory control model formulated in servomechanism theory 
by Vassian [6] and extended by Reilly [5] and Zehna [7], The 
techniques of linear system analysis as developed by Howard in 
[2] and the exponential smoothing/forecasting techniques of 
Brown [l] are heavily relied upon. The model as it is described 
in [7] will henceforth be referred to as the "basic model." 
Briefly, in the basic model we have demand as the input and 
inventory level, measured relative to some safety level, as the 
output. Inventory is replenished using a reorder rule which is 
based on a linear combination of past demand and inventory level 
feedback. Units of inventory are assumed to be issued instan- 
taneously and a fixed lead time is assumed for receipt of reple- 
nishment orders. In discrete variable servo theory notation the 
two equations which describe the basic model are 



(1-1) I-I.+e.-X and 

t t-1 t-l-T t 




where 1^ m on hand inventory relative to some safety level 
at time t. 



9 = reorder quantity computed at time t. 

t 



{a ) = a set of constants which weight past demands 
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{b^ } = a set of constants which weight past inventory 



levels. 

00 



Ic r ■> 

Let D(z) 5 ^ d k z • he the ^-transform of the sequence (d^j, 
k«o 



Then, the above equations become, in transform notation, 



$ . 



T+l 



(1-3) I(z) ■ 2I(z) + 8 9(z) - X(z) and 



(1-4) 0(z) - A(z)X(z) +B(z)I(z). 

It Is easily shown that 

,'T+l 

(1-5) I(z) - ^ ^ 

tfc r. ,-T+l 

i*4 - 4 b(z) 



X(t) 



Let 



d-6) S(z) “ 



. iT+1 

l4 A(z) -1 ] 

s T+l 
1-8 “ft B(z) 



be called the transfer function. Then. (1-5) may ba written 
I(z) - S(z)X(z) . 

Similarly, 

(1-7) 0(z) - [A(z) + B(z) S(z)] X(z) 



- s'(z) X(z) 

where S / (z) ■ Ca(z) + B(z) S(*)] 

the basic model is summarized in flow graph form in Figure 1, 

As Vassian suggests [6], we will let 
. -( 1 - 2 ) 



(1*8) B(z) 
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FIGURE 1 

flow graph diagram of the basic model 



It follows Chen that S(z) will have no finite poles. In that 
case, the denominator of S(z) will be finite for all finite values 

of s. 

> 

The following additional relationships, developed in [7], 
are listed for ready reference: 



(1-9) E(z) s eCi(z)] - S(z) m(z) 
where 

m(z) ■ E[X(z)] 

T+l (l-z T+1 ) 

(1-10) A*(z) - z A(z) A * X(z) 

1-z 



( 1- 11) A* ** i + \ X 

t+l-HT t+l+T L> t+l4T-j . We require that 

j“o 



t+l-fT 

(1-12) EiA* ] - Y m(j)->0 as 

t-J-l-KC u ‘ 

j~t*i'l 



When expected demand does not vary with time (E[X t ] « a) we have 
the following additional relationships* 
t 

/a V 7 k 

(1-13) a « cv ) 0 X , using exponential smoothing, 

t ^ t-k 

k=o 



/ or X(z) 

(l-H) a(z) = — — 

1 - P z 



(1-15) A* ~ (T-l-1) a , a forecast of demand during lead time 

t!l!T t 

plus one unit of time. 
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(1-16) A*(z) 



(1-17) A(z) = 

(1-18) S(z) = 

(1-19) I(z) = 

(1-20) 9(z) = 

Inventing the 

( 1 - 21 ) 9 = a 

t 

and 

(1-22) 9 = cr 

t 

Inverting the 

(1-23) I = a 
t 



T+l 

<* (T+l) z X(z) 
1 - 8 z 



<y (t+i) ( l-z) 

T+l 

(l-z )(l- 8 4 ) 



T+l 

or (T+l) 2 



1 - 8 k 

k 



T+l 

a (T+l) z X(z) 



1- 8 z 



a -jTS 

1-4 

& \ 



T+l) 

(1 z X(z) 



l-z. 



s'(z) X(z) 



« (T+l) (1 - 2 ) , ]> 
l-z' 



X(z) 



transform in (1-20), we obtain 



t t-1 

(T+L) } 8 k X - of (T+L) / 8 X , + X 

L, t-k t-l-k t 

k=o k«o 



(T+l) [X - a ] + X . 
t t-1 t 



transform in (1-19) yields 



t-T-1 

(T+l) l 8 k 
k=o 



X 

t -T-l-k 



T 




k=o 
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Computing the expected value of 1^ we obtain 

t-T-1 T 

(1-24) B "O' (T+l) £ e k m(t) - ^ 

k=o k-o 

t-T-1 

•r~t k 

“ or (T+l) a\ P • (T+l) a 
k=o 

r -ii 

- a (T+l) a - (T+l) a or, 

0 * 1 

t-T 

(1-25) E t - -a (T+l) & 

Thus, we have described the basic model and displayed the 
resulting mathematical relationships. In addition, for the 
case of constant mean demand we have derived explicit formulae 
for the expected on hand inventory and reorder quantities at 
time t when exponential smoothing is used to estimate the mean. 
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2. DIGITAL COMPUTER SIMULATION 



The simulation presented here was programmed in FORTRAN 
60 and run on the CDC 1604 serial number 1. The programs used 
are included in Appendix II in the order of discussion below. 

First, a deterministic simulation of the basic model with 
constant mean was made based on the relation 



letting a = .2, P - .8, and T - 4. in the first run m(t)-16 
and the results of 995 iterations are shown in Appendix I-A. 
It is seen, that to six decimal place accuracy, the expected 
inventory level reaches zero after 83 time periods. This is 
merely a simple verification of the analytical results shown 



In the second deterministic run, we let m(t)«t in the basic 
relation. Since Reilly [5] has shwn that first degree expo- 
nential smoothing yields biased results when the demand is 
linear we would expect this bias to be reflected in the expected 
inventory. In the case at hand we have 



t-T-1 



T 




t-T-1 



T 




Letting or “ .2, P » .8, m(t) ■ t and T « 4 



.i 



t-5 



4 




k=o 



k=o 
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t-5 t-5 t-5 v * 4 



^(V-sr^-r^ -svr* 

Iqssq k=o 



4 ^ 

k=o 



k«o k=o 



Izi^L 



t-4 



.2 j 



+ .8 



t-4 

■kl-J) 

.2 



t-4 t-5 

•(1-.8 ) - . 2(t-4) .8 



.04 



•5t + 10 



= -5t(.8) t-4 + 45(.8) t ' 4 - 4t(.8) t_5 + 16(.8) t ' 5 -35. 



Hence, lim E * “35 
t-? ® 



We can see from the results shown in Appendix I-B that in 
fact * “35 after 78 iterations again to six decimal place 
accuracy. So we have verified the analytical conclusions in 
[5] and [7], i.e., that the model is asymptotically unbiased 
(no steady state error in the expected inventory level) when- 
ever asymptotically unbiased forecasters are used end asympto- 
tically biased when correspondingly biased forecasters arc used. 
However, as a periodic review inventory model, the basic model 
may require an impractical number of periods to converge to an 
acceptable level. Conceding that convergence can be accelerate..* 
by adjusting a and T, we must face the fact that asymptotic 
unbiasedness may not be a practical measure of effectiveness for 
an inventory model. 
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The second simulation to be considered is stochastic in 



nature. Assume that demand is Poisson, i.e., 



P[X t =k] 



-x t k 



k! 



where X. ® X. . 



For ease of simulation a normal random number generator was 
employed taking advantage of the fact that for large X [ 3 ] 



-X 



e 



k! 




k-X + .5 

fT ay. 



k-X - .5 

rr 



The 1000 random numbers for the first two runs are listed in 

Appendix 1I1-A, while those used in runs three and four are 

shown in Appendices III-B and III-C respectively. 

In the first run X ■ 16 in the basic relation 
, t 

t-T-1 T 

i t -«(m>£ -I\-k • 

k=o k=o 



I thru I are listed from right to left in Appendix I-C. 

6 1000 

It is seen that there is some tendency toward stability through 
I ^ but after that time period we observe the effects of the 
random demand on the inventory level. It is now clear that 
carrying expected values through the model does not adequately 
describe the true situation. 

In the second stochastic, run, all parameters and variables 
were unchanged from the first run except that T was increased 
to T ■ 10. We see the resulting inventory values in Appendix I-D. 
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As expected, some increase in the variation of 1^ is observed, 
i.e., the variance of the random variable I fc is a function of 
the lead time of the basic modal. However, we also observe 
that is relatively insensitive to T in the sense that the 
order of magnitude of the change in T is not carried over to 
I t to an appreciable degree. 

In the third stochastic run, all parameters and variables 
were the same as in the first run except that X ■ 100. Natur- 
ally, to preserve our Normal approximationj the Poisson distri- 
bution variance was also increased to 100. Note that the random 
numbers generated and listed in Appendix IXI-B reflect this 
change. Now we observe the resulting values of 1^ shown in 
Appendix I-E. Analogously we see that the variance of 1^ 
depends on the variance of but is relatively insensitive. 

The order of magnitude of the increase in X is not reflected 

in I . 
t 

In the final run X t was allowed to vary with time, speci- 
fically X ■ 16 + t. From the analytical development we are 
led to expect bias since our forecast is biased for non-con- 
stant mean demand. In Appendix III-C observe the new set of 
random numbers generated with mean X = 16 + t and variance 
X^ ** 16 + t. The expected bias is apparent from an inspection 
of the Appendix I-F. Most values of I fc are negative. Again 
observe the relative insensitivity of I^_ to increasing the 
variance of the demand. 

In summary, the computer simulations have confirmed the 
theoretical results as expected./ However, in doing so, some 
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limitations of the basic model have been suggested. Asymptotic 
unbiasedness as a criterion appears to bo weak and probably 
impractical even when expected values of the random variables 
X and l t are considered. When {x..} is a sequence of realizations 
of a random variable from a probability distribution, the sequence 
£l t ) certainly does not appear to approach zero, Vassian's [6] 
apparent claim to the contrary. Zehna [7] clarified the theory 
in this area but the results of the simulation dramatically 
emphasize the inadequacy of the basic model when only expectations 
are considered. The choice of a safety level, or other managerial 
decisions, cannot be made on the basis of e[‘I J alone but Var[l t ] 
must be accounted for also. 
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3. DOUBLE ECHELON MODELS 



It is possible to formulate multiple echelon models in 
servo theory along the lines suggested by Howard [2]. It is 
desired to test the basic model of the previous sections in 
echelon. First, we construct a double echelon model using 
three basic models, two at the retail level and one at the 
wholesale supporting stock point position. Assume that con- 
sumer demand with constant mean values a^ and a^ is the system 
input at the retail level. In addition to the inventory level 
outputs at the retail stock points we have the reorder quantities 
6^(t) and 9^(t) serving as demand inputs to the wholesale stock 
point. There is a lag time which is constant for each retail 
stock point of and T^ respectively between the time of issue 
by the wholesale stock point and receipt by the retail points. 

Assume that the wholesale stock point has a sufficiently high 
safety level so that it is always able to fill orders from the 
retail activities. There are to be no interactions between the 
retail stock points. All three stock points are to use a single 
exponential smoothing forecaster which is of course asymptotically 

unbiased for constant mean demand. The lead time, T , for the 

3 

wholesale stock point is the delivery time from the manufacturer . 
Assume that all three stock points use the same value of the 
smoothing constant or. Hence this double echelon model describes; 
in a simple way, a common situation involving two echelons of 
inventory. The use of a common inventory control doctrine and 
the assumption of no interactions at the retail level are considered 
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